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We study stability and distortions of liquid crystal nematic order in a cell with a random hetero- 
geneous substrate. Modeling this system as a bulk xy model with quenched disorder confined to a 
surface, we find that nematic order is marginally unstable to such surface pinning. We compute the 
length scale beyond which nematic distortions become large, and calculate orientational correlation 
functions using the functional renormalization group and matching methods, finding universal log- 
arithmic and double-logarithmic distortions in two and three dimensions, respectively. We extend 
these results to a finite-thickness liquid crystal cell with a second homogeneous substrate, detailing 
crossovers as a function of random pinning strength and cell thickness. We conclude with analy- 
sis of experimental signatures of these distortions in a conventional crossed-polarizer-analyzer light 
microscopy. 

PACS numbers: 61.30.Dk, 61.30.Hn, 64.60.ae, 79.60.Ht 
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I. INTRODUCTION 
A. Motivation and background 

Over the past several decades, there has been con- 
siderable progress in understanding the phenomenology 
of ordered condensed states subject to random hetero- 
geneities, generically present in real materials 0, 0. 
These include "dirty" charge-density waves 0, super- 
conductors [J, and magnets, as well as superfluid helium 
[j| and liquid crystals p-Q m the random environment 
of aerogel and other porous matrices. 

Much of the detailed understanding came from the 
analysis of the xy-, O(N)-, and related smectic @, 0] 
random-field models, in the pioneering works by Larkin 
E3] i Fisher and Nattermann [12[ , and extensive 
subsequent studies by Le Doussal and co-workers [ilj - fl6| 
using a combination of replica variational and functional 
renormalization group (RG) methods. At low tempera- 
ture and weak disorder (neglecting enigmatic effects of 
topological defects that may become important on much 
longer scales), the state is characterized by elastic distor- 
tions with universal power-law correlations controlled by 
a nontrivial zero-temperature fixed point, the so-called 
"Bragg" (elastic) glass 0, QtHU . 

With the exceptions of the pioneering surface disorder 
study by Feldman and Vinokur [21[ and our recent exten- 
sion to liquid crystal cells [22J , all of the theoretical foci 
have been on the bulk heterogeneity, where the disorder 
extends over the full sample. However, there are many 
realizations in which, instead, random pinning is confined 
to a subspace, e.g., a surface of the sample. A techno- 
logically relevant example illustrated in Fig. [T] is that of 
a liquid crystal cell (e.g., of a laptop display), where a 
dirty substrate imposes random pinning, that competes 
with liquid crystal ordering. In this paper, as a signif- 
icant elaboration of our earlier work [221 ] . we present a 
detailed analysis of such a cell. 

A commonly observed Schlieren texture [23] is a man- 
ifestation of such surface pinning in nematic cells. Re- 
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FIG. 1: (Color online) Nematic (top) and smectic (bottom) 
liquid crystal cells of thickness w with both substrates favor- 
ing a planar alignment. The front substrate favors a random 
heterogeneous in-plane orientation, while the back substrate 
is rubbed and therefore picks out an ordering orientational 
axis. 



cent studies also include photo- alignment and dynam- 
ics in self-assembled liquid crystalline monolayers |24| , as 
well as memory effects and multistability in alignment of 
nematic cells with heterogeneous random anchoring sub- 
strate [25[. The existence of the corresponding phenom- 
ena in smectic liquid crystals has been recently revealed 
in ferroelectric smectic-C cells in a book-shelf geometry 
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(26| . This latter system exhibits a long-scale smectic lay- 
ers distortion, demonstrated to be driven by collective 
random surface pinning, and awaits a detailed descrip- 
tion. Although here we do not provide an analysis specific 
to these experiments, we expect that results presented in 
this paper are a first step toward this goal. 

We conclude the introduction below by summarizing 
our main results. The rest of the paper is organized 
as follows. In Sec. [n] as a simplest description of ne- 
matic liquid crystal cell, we present a model for an xy 
system ordered in a homogeneous bulk, but randomly 
pinned on a surface. In Sec. IIIII we analyze this model 
within the random-torque Larkin approximation, valid 
on short length scales. To study the physics on longer 
length scales we employ the functional renormalization 
group (FRG) in Sec. IIV A[ combining it with matching 
methods in Sec. IIV Bl to compute the asymptotics of ori- 
entational correlation functions. We briefly consider the 
limit of strong surface pinning in Sec. |Vj In Sec. I VII 
we discuss application and extension of our findings to 
nematic and smectic liquid crystals, and analyze exper- 
imental signatures of heterogeneous surface-driven dis- 
tortions in a crossed-polarizer-analyzer light microscopy. 
We conclude in Sec. IVIII 



B. Summary of the results 

In this paper, we study stability and distortions of a 
liquid crystal order in a cell with a random heterogeneous 
substrate. We model such and related systems as a bulk 
d-dimensional xy model with quenched random-field dis- 
order confined to a (d — l)-dimensional surface 27|, and 
analyze its stability to surface pinning and compute cor- 
responding correlation functions. As we will show, many 
of the tools developed in the context of the bulk disor- 
der can be taken over to this ubiquitous surface-pinning 
problem. 

Since the random pinning is confined to a surface, one 
might a priori expect its effects to be vanishingly weak 
compared to the ordering tendency of the homogeneous 
bulk, and thus the xy order to be stable to weak surface 
disorder in any dimension. Our finding contrasts sharply 
with this naive intuition [28j . Namely, our key qualitative 
observation is that the xy order in such d-dimensional 
system with d < di c , where 



dir. — 3, 



(1) 



is always destabilized by arbitrary weak random surface 
pinning [2ll l22j|. For d > 3, the xy order requires a finite 
threshold of heterogeneity to be destroyed. Thus, while 
the lower-critical dimension of di c = 3 for the system 
with surface disorder is indeed reduced relative to that 
of the bulk pinning problem, (characterized by = 4 

[llllilj]); surface heterogeneity has a qualitatively strong 
effect on the bulk ordering even in the thermodynamic 
limit. 



Above observation can be simply understood from 
a generalization of the Imry-Ma argument [29[ to the 
surface-pinning problem [2TL [22| ■ For an ordered region 
of size L, the interaction with the surface random field 



can lower the energy by an order of E p i n ~ V p yN p ~ 

A l/2 i( rf-l)/2 ; where y p 

is a typical pinning strength with 

zero mean and variance Af ss V 2 /^ -1 (£0 is the pinning 
correlation length) and N p is the number of surface pin- 
ning sites. Since a surface distortion on scale L extends 
a distance L into the bulk, the corresponding elastic en- 
ergy cost scales as E e ~ KL d ~ 2 , where K is the elastic 
stiffness. By comparing these energies, it is clear that for 
d < 3, on sufficiently long scales, L > £l ~ (K 2 /Af)^^, 
the surface heterogeneity always dominates over the elas- 
tic energy, and thus on these long scales always destroys 
long-range xy order for an arbitrary weak surface pin- 
ning. 

A more detailed analysis extends the argument to 
three dimensions (3D). The corresponding surface Larkin 
length scale [l(|, beyond which the T = orientational 
order on the random surface (z = 0) is destroyed, is given 
by 



for d = 3, 



(3-d){2ir) 2 K 2 



Cd-xAf 



for d < 3, 



(2) 
(3) 



with a is a microscopic cutoff of order of a few nanometers 
in the context of liquid crystals, set by the molecular 
size, c = 8tt 3 , and C d _i = 2 2 - < V 1 - d V 2 /r(^i) as in 
Sec. HVA1 

At a distance z»a away from the substrate (into the 
homogeneous bulk), the orientational order distortions 
(characterized by mean-squared fluctuations of 4>) across 
a region of size £l decay according to 



1 - If (lnfe + l-7), 2z«£l 



£L e -2z/£ £ 



for d 



l-T{d-2)(^f- d , 2z«£ L 



(3-d) 



2z c ' 



2z> i L 



for 2 < d < 3, 



(- , _ ln(2z/q) „ f 

- t /2 n i ?i/a V/, ford = 3, (4) 

I lnfa/a) e > lZ » 5£ 

where 7 ss 0.58 is the Euler's constant. Thus we find 
that the orientational order distortions induced by the 
heterogeneous surface penetrate a distance £l into the 
bulk. 

The physics is more complex on in-plane (x) length 
scales longer than £l, where the random-torque Larkin 
approximation is invalid and the <f> nonlinearities of the 
pinning potential, V(0,x), must be taken into account. 
To access these longer scales, we employed FRG and 
matching methods. On the substrate (z = 0), for d < 3 
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we find (q is an in-plane wave vector) 



for d < 3, 



(5) 



where c(d) is a universal number [given in Eq. (I89|) ]. that 
in the physically relevant case of d = 2 is given by c(2) = 
At the lower-critical dimension, d = 3, the distortion 
variance is given by 



q » Vfc 



A/ 1 

~ ~l~ q2 ln(ga) ' 9<!/Ci 



for d = 3. (6) 
These give the correlation function C(x, 21,22) = 



((</>(x, zi) — </>(0, z 2 )) 2 ), where (• • • ) is a thermal average 
and 777 is a quenched disorder realization average, with 
latter dominating the former at low temperatures. We 
report its full spatial dependence in Sees. Mil and IIVBI 
and in Figs. [10] and [T3l On the substrate (z — 0), its 
asymptotics in two dimensions is given by 



0,0) 



for with finite z, the correlation function is given 

by 



C&\x,z,z) « 6 2e - 2 ^+^ln 



1 



-(8) 



in which 6 2 is a weak function of 2z . In the physically 
most relevant three dimensions, on the substrate, it is 
given byHJll 



C^(x,0,0) 



Sir 2 



2si In 
9 



ln(x /a) 
ln(Wa) 



( 9 ) 

At large x 3> £l with finite z, the correlation is given as 



2_7T 2 

9 



16z a In (2z/a) ' 



< ac < 2z, 
(10) 



in which 63 is also a weak function of 2z/£l- The ex- 
ponentially decaying parts of both two-dimensional (2D) 
and 3D results are contribution from short scales below 
the Larkin length, for z 3> £z,. 

With an eye to liquid crystal cell applications, we also 
analyzed a nematic cell of finite thickness w, with (as 
above) a heterogeneous bottom substrate, and a top sub- 
strate with a homogeneous Dirichlet or a homogeneous 
Neumann boundary conditions. Not surprisingly, this 
presents new crossover as a function of the ratio, w/^l, of 
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FIG. 2: (Color online) Orientational order parameter tl>(w, 0) 
(controlling light transmission through a liquid crystal cell) at 
the random pinning substrate of a 3D Dirichlet cell of thick- 
ness w. 



the cell thickness w to the Larkin length in a (infinitely) 
thick cell. 

For the Larkin length in the Dirichlet cell, illustrated 
in Fig. [5] we find 



Si 



£,L, £l < w 



(11) 



where = a^w is the crossover "bulk" Larkin length 
beyond which effects of the finite cell thickness become 
important, c 2 — l,a 2 ss 1.71, 1/3 = 1 and C3 ~ 0.79, 03 
1.23,^ = 1/2. 

Above behavior in the Dirichlet cell is a manifestation 
of a crossover as a function of cell thickness (or equiva- 
lently disorder strength ) from a weakly ordered state for 
a thick cell (and strong disorder) to a strongly ordered 
state for a thin cell (and weak disorder). The crossover 
is more clearly reflected in the surface xy orientational 

that in d < 3 is given by 



order parameter ip = (e 1 ^ 

^ - a (w/£ L f- d 
1Pd<3 ~ < ..ft. \ r >i 



thin cell, w -C £l 



thick cell, w» a, 



(12) 



and in 3D the orientational order parameter, illustrated 
in Fig. [21 is given by 

^3D ~ S „ 



ln(q/a) 
ln(iti / a) 



thin cell, w <C £l 

V3D (13) 

, thick cell, w £l, 



where i]^ = (3 — d)7r 2 /18, 773.D = 7r 2 /18 are universal 
exponents [given in Eqs. (|M| and (llOlj) ] and a = 2tt 2 , 
T]l = 2ir 2 / ln(£i/a) are nonuniversal constants. 

In contrast, for a cell with a Neumann boundary con- 
dition on the top substrate, the Larkin length, illustrated 
in Fig. [51 is given by 



AW 



w 1d 



y/2\n(1.2^ L /w), 



3-d 
5-d 



d<3 



w > £l 



(14) 
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with = and the state is disordered, i.e., tjj = for 
arbitrary weak random pinning. 

We expect these crossovers as a function of to 
be relevant to understanding the ordering in real liquid 
crystal cells. They should be accessible experimentally 
in a setup of the type used in Ref. [Hj]. 



II. THE MODEL 



As a "toy" model of a nematic liquid crystal cell with a 
dirty substrate and thickness w, illustrated in Fig. [3l we 
employ a d-dimensional surface random-field xy model 
characterized by a Hamiltonian 



H= I d 

' — oo 



d-1 



c I dz 

o 



K 



(V0(r)) 2 - V[0(r),x]5(z) 



(15) 



In above 0(r) is the xy-field distortion at a point r = 
(x, z), the random pinning potential V[0(x, z), x]<5(z) 
with a 27r periodicity of 0, characterized by zero mean 
and Gaussian distribution with a variance 



V((f>, x) V(0', x') = i?(0 - <j>')8 d - 1 (j 



(16) 



is confined to the bottom substrate at z = 0, and we im- 
pose either a Dirichlet [0(x, z)\ z — w = 0] or a Neumann 
[d z (j)(x., z)\ z=w — 0] boundary condition on the top ho- 
mogeneous substrate at z = w . 

A more realistic model must of course include nonpla- 
nar director distortions, characterized by an additional 
polar angle, as well as point and line (strength 1/2 discli- 
nations) topological defects allowed by the 3D headless 
nematic director field, n. We discuss the effects of these 
additional ingredients in Sec. I VII The long-scale behav- 
ior of the coarse-grained periodic (period 2tt) variance 
function R(<fi) characterizes low-temperature properties 
of our system and will therefore be our main focus. 



A. Dimensional reduction 

Because the random pinning potential in the Hamilto- 
nian, (fl~5j) . is confined to the bottom substrate at z = 0, 
no nonlinearities appear in the bulk (0 < z < w) of 
the cell. Consequently, as in other similar problems 
[30l l3l| . it is possible and convenient to focus on the 
random substrate and exactly eliminate the bulk degree 
of freedom of 0(x, z) in favor of the random substrate 
field </>o (x) = 0(x, z = 0). This can be done via a con- 
strained path-integral by integrating out 0(x, z) with a 
constraint 0(x, z = 0) = 0o(x), thereby obtaining an ef- 
fective (d — lVdimensional Hamiltonian for 0o(x) (30j . 
Equivalently [HI H|[ , we can eliminate 0(x, z) by solv- 




FIG. 3: Schematic illustration of the orientational distortions 
in an xy- model system confined to a thickness w cell with ran- 
dom pinning on the bottom substrate and Dirichlet boundary 
condition (rubbed substrate) on the top substrate. A Larkin 
domain (defined by the faint dashed curve) of size £l is also 
depicted. 



whose solutions for the boundary conditions of interest 
are obtained by elementary methods 



<^)(q,z) = O (q) 



Pl(n = tn] smh[ q (w-z)) 

sinh (qw) 

.(AO/ \ , , ,cosh[g(w - z)\ 
<p^'(q,z) = 0o (q)- 



w — y 00, (18) 
Dirichlet, (19) 

Neumann. (20) 



cosh (qw) 

We summarize above three boundary condition cases by 

(a) (q,z) = o (qV (a) (9^), (21) 

where mode functions ^^(q, z) are implicitly defined by 
Eqs. 

Substituting these into the Hamiltonian, Eq. (|15p . we 
obtain the (d — l)-dimensional (surface) Hamiltonian 



H, = 



''" \^lt>|0 o (q)| 2 - / d d - l xV[Mx),*}, 



(27T) 



(22) 



which characterizes the behavior of the field 0o(x) con- 
fined to the random substrate at z — 0. In the equation 
above, label a ranges over "free" (w —> 00), Dirichlet, and 
Neumann boundary conditions on the z = w substrate, 
with corresponding kernels given by 

= Kq, w^oo, (23) 

= Kqcoth(qw), Dirichlet, (24) 
r^> = Kqt&nh(qw), Neumann. (25) 



ing the bulk Euler-Langrange equation V 0(r) — 0. To . , „ . ,, . , ^. . , , , mCDK , 

.,? , „ ■ , r ,/ \ i , ■ • As expected, the hmtc-thickncss Dirichlet 1; ) and 



this end, we Fourier transform 
an equation for 0(q, z) = J d d ~ 1 



x, z ) over x, obtaining 
x<p(x, z)e _lq ' x , 



d 2 0(q,z)- 9 2 0(q,z) =0, 



(17) 



Neumann (Tq^^) kernels for w-ioo reduce to the case 
of the free kernel, Kq. The q nonanalyticity and long- 
wavelength stiffening (relative to the bulk Kq 2 kernel) of 



5 



the latter arises due to a mediation of surface distortions 
by long-range deformations in the bulk of the cell. In the 
opposite limit of a thin cell and long scales, as expected, 
the Dirichlet kernel reduces to a "massive" one K/w, 
and the Neumann kernel simplifies to Kwq 2 of an ordi- 
nary surface (without a contact with the bulk) xy model. 
The advantage of the dimensional reduction above is that 
formally the formulation of the random surface problem 
becomes identical to that of the extensively studied bulk 
random pinning [TTl-[l4| but in one lower dimension and 
with a modified long-range elasticity. 



B. Replicated model 

To compute self-averaging quantities, (e.g., the 
disorder-averaged free energy), it is convenient (but not 
necessary) to employ the replica "trick" [34], which al- 
lows us to work with a translationally invariant field the- 
ory at the expense of introducing n replica fields (with 
the n — > limit to be taken at the end of the calculation). 
For the free energy, this procedure relies on the identity 
for the ln(x) function 



F = -Tin Z = -T lim 



Z" 



1 



n->0 n 



(26) 



After replicating and integrating over the random poten- 
tial V[<^,x] using Eq. (fT6l) . we obtain 



Z" 



-m 



t>S]/T 



(27) 



The effective translationally invariant replicated Hamil- 
tonian [</>q ] is given by 



1 71 r 

-^fE/ ^-^[^(x)-^(x)]. (28) 



We will use this Hamiltonian, (|28p. in our subsequent RG 
analysis of the system. 



III. LARKIN ANALYSIS 

A. Random torque model 

As with the bulk quenched disorder, the nontrivial na- 
ture of the surface-pinning problem is encoded in the 
nonlinearity in <fio{x) of the random surface potential 
V[cf>o(x),x] in H s , Eq. ([22]) . However, in an approxi- 
mation first employed by Larkin Q (that now bares his 
name), for small <f>o distortions we can Taylor-expand the 



random potential to linear order (35) in 0o 

H * L) * I (|r)^5 r * a) l<Mq)| 2 - J rf d ~Wx)<Mx) 

* / [^(°»ko(q)| 2 -r^q)^(q) 



(29) 



obtaining a harmonic Hamiltonian characterized by (j)Q- 
independent random surface torque 



t(x) = 90 o V[0 o (x),x] 

= no, 4 



1=0 



(30) 



The random-torque inherits Gaussian statistics from that 
of the random potential V[(/>o(x),x], with its variance 
given by 



r(x)r(x') = -i?"(0)^ 1 (x-x') (31a) 
= A/^-^x-x'). (31b) 



B. Correlation functions 

The resulting surface Larkin model (valid on scales 
shorter than £l) is quadratic in <fio(x) and can therefore 
be analyzed exactly by standard methods. The basic 
quantity of primary interest is the Fourier transform of 
the disorder- and thermally averaged two-point correla- 
tion function 



Mq<> = (^q>(^') + ((0q-(0q))(^-(^») 

= [C A (q) + C T (q)}(27r) d - 1 d d - 1 (q + q'), (32) 



with 



A i 



r. 



(a) 



2 - 



a 



T, Larkin 



T 



r 



(o)' 



(33) 
(34) 



with the Larkin approximations for the T = distor- 
tions of </>o( x ) and its thermal fluctuations about this 
pinned ground state, respectively. Because of the respec- 
tive structures of Ca(<?) and Cx(q), at long scales (small 
q), clearly the quenched-disorder-driven distortions in the 
ground state dominate over small thermal fluctuations. 
Thus, for the remainder of the paper, we will focus on 
the system at T = 0. 

The real-space correlation function 



C(x,z,z')-((0(x,z)-0(O,z')) 2 ), 



(35) 



which describes surface and bulk distortions that fol- 
lows via a Fourier transform of Eq. (|33[) . combined with 
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Eqs. (|18p - ([20|) . Specializing for simplicity to the case 
z = z' , we compute 



Ci Q) (x,z,z) « 2A f J ^^(l-cosq.x) 



^ a \q,z) 



(36) 



for the three boundary conditions of interest, with the 
subscript "L" denoting that here the validity is limited 
to the Larkin approximation. These correlation functions 
are illustrated in Figs. HI [51 and [TO 

For the infinitely thick cell (w — > oo), the asymptotics 
(for a <C x <C £l) is given by (see Appendix IC 2\ 



C 



(°°) i \ 



2A 



(1 — cosq • x)e 2qz 



x, z = 

' 2 e~ 2z ^ L , z > a; 



for x < £ L , 



2z£ £ 

(t) 

3: 2 (2 Z ) 1 



3-ci 



for d 



z = 



for d < 3, 



In (x/ a) 
In (a /a)' 
(S) 2 (l+lf) - 



ln(?t/a) 



for d = 3, 



(37) 



and by definition vanishes as x — > 0. 

For a cell of thickness u> with a Dirichlet and Neumann 
boundary conditions on the top substrate, the correlation 
functions are given by 



Cf(x,v) 



<#°(w) 



2A 



^ 2 

2A, 



if 2 



(1 — cos q • x) sinh 2 q(w ~ z) 
q 2 cosh 2 qw 

(38) 

(1 — cos q • x) cosh 2 q(w — z) 



q 2 sinh 2 qw 



(39) 



For a thick cell (w > £l), these can be approximated by 
the infinite cell results given above. For a thin cell (w < 
£l) with a top Dirichlet substrate, their asymptotics is 
given by 



Cf } (x,z,z) 



(3-d)(fL)3-1, 



z <C w <C x 



( 3 ~ d ) Jd Z -lg3-d , £<W<X 



for d < 3. 



(40) 



The asymptotics for <i = 3 and for the Neumann bound- 
ary conditions are more involved and are best evaluated 
numerically. They are displayed in Fig. 01 



Ci a) (x,0,0) 




Q"(x,0,0) 



FIG. 4: (Color online) Numerically evaluated Larkin (a; < £l) 
correlation functions for a finite-thickness cell w, with the 
Dirichlet and Neumann top substrates in 2D (top) and 3D 
(bottom). In both figures, the thicker black curve corresponds 
to the infinite cell thickness, above it are the results for the 
Neumann substrate (with w/£l = 1,0.5,0.1, bottom to top), 
and below it are the results for the Dirichlet substrate (with 
w/^l = 1,0.5,0.1, top to bottom). As expected, the Dirichlet 
(Neumann) homogeneous substrate reduces (enhances) orien- 
tational distortions of <f>, with its influence growing with re- 
duced cell thickness and weaker disorder, controlled by 



C. Limits of validity: Larkin length 

An important feature of the Larkin approximation is 
that it can be used to estimate its own range of validity. 
The Larkin model (|29)) breaks down when the Taylor 
expansion in (f>o that led to it becomes invalid. This can 
be estimated by looking at the mean-squared distortions 
of 0o, given by 



A f 



,(a) 



(41) 



Focusing first on an infinitely thick cell, characterized by 
Tq°°^ = K q, we note that for d < di c = 3, mean-squared 
fluctuations of (x) diverge at long length scales. Thus, 
we find that plBif 

di c = 3 (42) 



7 



is the lower-critical dimension for the stability of the xy 
order in the presence of a random surface quenched pin- 
ning. This value of di c = 3 is to be contrasted with the 
^{buik) _ Q £ a g y S t em (spontaneously breaking 

continuous symmetry) subjected to bulk disorder. Thus, 
for d < 3, we can define the (so-called) Larkin length, £l 
as the scale at which O distortions grow to order 1 (that 
we take to be 2ir, for concreteness). That is, 



«p 2 (0,z)> 



^(x)>=(2tt) 2 



A/ 

K 2 



l ' a d d - l q 1 



V6 



(27r) d -! q 2 ' 



(43) 



which gives (see Appendix |A"|) the Eqs. © and ([3)) 

£ L = ae cK ^ A f, ford = 3, 

r (3 - 4(2^^1 A 



with a a microscopic cutoff of order of a few nanometers 
in the context of liquid crystals, set by the molecular 
size, c = 8tt 3 , and C d _i = 2 2 -<V 1 - d )/ 2 /r(^) as in 
Sec. HVA1 

Since the orientational order parameter 



P -W)/2 



(44) 



(somewhat crudely assuming Gaussian correlations in (f> 
in the second line above) decays with the growing 4>rms, 
the Larkin length, is the scale beyond which the ori- 
entational order falls off significantly. 

As discussed in the Sec. HI the mean-squared distor- 
tions of <?!>(x, z) are suppressed away from the random 
substrate, within Larkin approximation on scale £l, for 
w — y co decaying with z»a according to (see Appendix 

leu 



£l ^ m 2z 



li p -2z/£ 



2z 

2){f L ) 



l- 7 ), 2z«£l 
^ 3 - d , 2z«£ L 



(3-d) 



£ip-2z/? £ 
2z C 



1 - 



ln(2z/q) 



2z > Cl 

2z<a 



for d 



for 2 < d < 3, 



2z>a 



for d = 3. 



(45) 



The complete numerically-evaluated behavior is dis- 
played in Fig. [5] Thus, the orientational order heals on 
the scale of the Larkin length, into the bulk. 

In the interest of applications to liquid crystal cells, 
we generalize this analysis to a system with finite thick- 
ness, w (along z). For a homogeneous Dirichlet boundary 
condition on the top substrate of the cell, we have 



w) - w - R2 j (27r)d _ 1 



q tanh qw 



(46) 





2.0 



FIG. 5: (Color online) Mean-squared distortion within the 
Larkin approximation as a function of for 2D and 3D 

systems. 



which gives (see Appendix IA 3[) Eq. (TTTj) 

U « w 

with t;* L = adW the crossover "bulk" Larkin length, C2 = 
1,<22 ~ 1.71,1/a = 1; and C3 rj 0.79, 03 s» 1.23, t/3 = 
1/2. The complete behavior of C^ '^' m 2D and 3D is 
illustrated in Fig. [5] 

This behavior is quite consistent with physical expec- 
tations. For a cell thicker than the bulk Larkin length, 
£l (thick cell and/or strong disorder), there is little im- 
pact of the top substrate on the range of the finite xy 
order, dominated by the random lower substrate. How- 
ever, for a thin cell and/or weak disorder, such that 
the bulk Larkin length extends across the cell thickness, 
the Dirichlet substrate effectively enforces the xy-order 
alignment across the cell, suppressing </>g ms below 2-tt and 

(T>) 

thereby driving the cell Larkin scale, £ L , to diverge. 

It is important to emphasize that this divergence is 
not an indication of a sharp transition. One signature of 
this is the fact that the "crossover" scale is a func- 
tion of a relatively arbitrary constant [taken here to be 
(27r) 2 ] in the definition of the Larkin length. Rather, the 

divergence of £^ , Eq. (fTTj) . is a signal of a crossover 
from a weakly xy-ordered state (at strong disorder and 
thick cell) for (1 < w to a strongly icy-ordered state (at 
weak disorder and thin cell) for £l ^> w. In both lim- 
its, the aligning Dirichlet substrate dominates over the 
random one, leading to a long-range xy order. We will 
support this assertion with a detailed FRG calculation in 

Sec. EES 

For a cell with a Neumann boundary condition on the 
top substrate, the Larkin length, illustrated in Fig. [6] is 
defined by 



(AO 



(27T) 



K 2 



1/a d d- 



1 



/( m {2ir) d - 1 g 2 ta nh 2 



qw 
(47) 
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FIG. 6: (Color online) Larkin lengths with different bound- 
ary conditions on the top substrate, for d = 2 (top) and d = 3 
(bottom). For thick cells (w S> the asymptotic behavior 
£L D ' Af ' ~ is as expected. For longer £&, i.e., weaker surface 
pinning and thinner cell, f^ 13 ' i s always longer (shorter) 

than £l for the infinite-thickness cell, consistent with the or- 
dering (disordering) effect of the finite cell thickness. For the 
Dirichlet homogeneous substrate, the effective Larkin length 
diverges at £l/w « 1.71 for d = 2 and £,l/w « 1.23 for d = 3, 
with asymptotics quoted in the text. 



This gives 



AN) 



with 



w 



Id 



V21n(1.2&/«;), 
(fEi)^(^) 1 ^. 



d= 3 
d < 3 



Id = 



(48) 



(49) 



Strongly contrasting to the case of the Dirichlet top sub- 
strate, this result is again consistent with our expecta- 
tions. Namely, that for a thin cell, w <C £l, such Neu- 



mann cell becomes a thin (d— l)-dimensional film, reduc- 
ing to a (d — l)-dimensional xy system pinned by (d— 1)- 
dimensional, i.e., bulk disorder, characterized by the ex- 
pected lower critical dimension of 5 (that is, di c — 1 = 4). 
Thus on length scales exceeding the Larkin lengths, t^£\ 
<po distortions become large and lead to a breakdown of 
the Larkin model, (f2U)) and of the predictions of correla- 
tion functions [e.g., Eqs. (|3"T|) and (|4"0)> ] calculated from 
it. 



IV. PHYSICS BEYOND £ L 

On length scales longer than the crossover scale ^\ 
large 4>q distortions are in the nonlinear regime and the ef- 
fects of the random potential V[0o(x),x] must be treated 
nonperturbatively. As with the bulk disorder problems, 
this can be done systematically using a renormalization- 
group analysis [ill ITij in an expansion in e = di c — d = 
3 — d about the lower-critical dimension. 



A. Functional renormalization-group analysis 

It is convenient to work with the translationally invari- 
ant replicated Hamiltonian, Hs r \ Eq. (|2"8)) . We employ 
the standard momentum-shell RG transformation [36] by 
separating the xy field into long- and short-scale contri- 
butions according to <^q (x) = (x) + (x) and per- 

turbatively in nonlinearity R\4>q (x)— 0q ( x )] integrate out 
the high wave vector piece 0q>(x) that takes support in 
an infinitesimal shell A/6 < q < A = l/a, with b = e se . 
We follow this with a rescaling of lengths and the long- 
wavelength part of the field 



x 

q 

0o<(fcx') 

<h<(d/b) 



6x', 
b d-x 



(50a) 
(50b) 
(50c) 
(50d) 



so as to restore the UV cutoff back to A. In Eq. (I50c|) . we 
made a convenient choice of a zero scaling dimension for 
the real-space field 0o( x ). This is dictated by the conve- 
nience of keeping the periodicity of the disorder variance 
R{(j>) fixed at 2tt. 



1. Infinitely thick cell: w — > oo 

We first focus on an infinitely thick cell, defined by 
w >• £l- The above rescaling leads to zeroth order RG 
flows of the effective couplings [37[ that for a thick cell is 
given by 



K(b) 
R{4>,b) 



b d - 2 K , 
b d - x R{cj>), 



(51) 
(52) 
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indicating that in d > 2, the effective strengths of both 
elastic and pinning energies grow at long scales relative to 
the thermal energy, T. This is a reflection that in d > 2, 
the physics is controlled by the zero-temperature ground- 
state competition between elastic and pinning energies, 
at long scales both much larger than the thermal energy. 
Equivalently, to emphasize this physics, we can rescale T 
according to 



T(b) = b-V-^T 



(53) 



so as to keep the elastic energy fixed at order 1. With 
this convenient rescaling convention, the measure of the 
effective pinning strength grows according to 



R(4>,b) = b 3 - d R((j>), 



(54) 



modified by a factor (T(b)/T) 2 = b 2( - d - 2 ~> relative to that 
in Eq. ([55) due to the factor of 1/T 2 in iji r) /T, Eq. {25]). 
Equivalently, without the rescaling of T, the dimension- 
less combination that arises in the coarse-graining analy- 
sis is given by R(4>) /K 2 , and its zeroth order flow is given 
by Eq. dMJ). 

In either convention, we find that for d < 3, the influ- 
ence of the random surface pinning grows at long scales 
relative to the elastic energy, consistent with the scaling 
and Larkin analysis that gave di c = 3. 

The statistical symmetry[38j of the bulk Hamiltonian, 
H (fTB"]) . under an arbitrary local rotation <fr(r) — > </>(r) + 
X(r) guarantees that the flow of K (b), Eq. ([ST]) , and 
equivalently, the thermal exponent 



9 = d- 2, 



(55) 



are exact, i.e., do not experience any coarse-graining cor- 
rections. This can equivalently be seen from the repli- 
cated Hamiltonian (|28[) . where the pinning nonlinearity, 
R[(f>o (x) — 0q (x)] depends only on the difference between 
different replica fields, i.e., independent of the "center of 
mass" field X)q=i^o- That is, the only nonlinearity in 

(r) 

Hs exhibits a symmetry of a replica-independent local 
rotation (r) — > (j)^ (r) + x(r) and under coarse graining 
can therefore only generate terms that also exhibit this 
symmetry. Thus, it cannot generate a correction to the 
elastic term that clearly lacks this symmetry, implying 
that K is not renormalized by the pinning disorder. 

An important consequence of the periodic nonlinearity 
R(<p) and the effective zero-temperature physics, first em- 
phasized by Fisher (llj . is that all monomials or (equiv- 
alently) harmonics in the expansion of R(4>) are equally 
relevant in d < di c . Thus, a functional RG analysis that 
follows the coarse-graining flow of the whole function 
R{4>) is necessary. The method is by now quite stan- 
dard [HI EH-dl] and is straightforwardly adapted to the 
surface-pinning problem, characterized by Hs , Eq. (|28|) . 

We limit the FRG analysis to one-loop order, perform- 
ing the momentum-shell integration over the high-wave- 
vector components 4>o> perturbatively in the nonlinearity 



R[<Pq (x) — </>q (x)]. We find that the change in the Hamil- 
tonian due to this coarse graining is given by 



a <+0o>])>-- - , 



(56) 

where H p [(/)q] is the nonlinear pinning part of the Hamil- 
tonian H { s r \ Eq. pg]). 

n . 

ff P = -^E dt^xRffiW-tfiW], (57) 

and the averages over short scale fields, <^o>, above are 

performed with the quadratic (elastic K) part of H^ r \ 
The superscript c denotes a cumulant average, (H 2 ) c = 
(H 2 )-(H p ) 2 . 

To lowest order in R(cf>) (dropping a constant term) we 
find that SH^lcft^] is given by 

a,P x 

« -j^E / R"^<-4<)m>-^>?)> 

a,P x 

* -^r^E/ (58) 



'1 L 1 a, 

which when compared to the definition of H p gives 
6RW(<p)K5lC d - 1 A d - 2 ^R"(4 > ). 



(59) 



In above, the prime indicates a partial derivative with 
respect to <f>, and C d = S d /(2n) d = l/[r(d/2)2 rf ~ V*/ 2 ], 
with Sd the surface area of a ci-dimensional unit sphere. 
The contribution to second order in R{<j>) is given by 



6 s2 8T 3 2 



J_l y 



i?"[^( Xl )-^( Xl )] 

xi?"[0-(x 2 ) - ^(x 2 )]/^(x 1)X2 ) , (60) 



where 



d 2 = ^((^(^-^(xi)) 2 ^^)-^.^)) 3 )". 

(61) 

Keeping only the most relevant (two-replica) terms and 
comparing to H p , we obtain 



SR<- 2 \<t>) a S£g 2 ( l -R"{4>)R"{4>) - R"(4>)R"(0)) , 



where the constant 172 is defined by 



(62) 



6£g 2 



1 



K 2 



(63) 
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Combining the first and second-order contributions to 
R(4>), Eqs. and (|()2"j). with the length and field rescal- 
ings, Eqs. (fSTj) and (|52|) . we obtain the FRG flow equation 

diR(4>) = eR{<j>) + C d - l A d - 2 ^R"{<i>) 

+lR"(<P)R"(<t>)-R"(cf>)R"(0), (64) 



where 



(65) 



is the dimensionless measure of surface disorder. 

(i) 2 < d < 3. Because as noted above T/K flows to 
zero as b 2 ~ d — 6~ e , for d > 2 (independent of the rescal- 
ing convention), i.e., the system is described by the zero- 
temperature fixed point, the second term on the right- 
hand side in Eq. (|64|) can be neg lected near d = 3, and 
FRG equation reduces to [U [jj, [H [H| 

d t R{<t>) = eR(cP) + ±R"(cf>)R"(<f>) - R" (cf>)R"(0). (66) 

We note that aside from constant prefactors in the defi- 
nition of R(4>, £) and the reduced lower-critical dimension 
giving e — 3 — d, the flow equation for the dimensionless 
disorder measure R(4> : £) is identical to that of the bulk 
pinning problem [111 Il4j . Consequently, the long-scale 
properties of the low-temperature phase are described by 
the same fixed point function, 



R'M 



7T" 

18 



(67) 



periodically extended (period 27r), with the minimum at 
the cusp given by R"(0) = — ^f-- 

(ii) d — 3. Temperature remains irrelevant at d = 3 
(as for any d > 2), allowing us to continue to work at 
T = 0. At this lower-critical dimension, e = 3 — d = 
and the flow equation reduces to 

d e R(cj>) = ±R"(4>)R"(4>) ~ i?»i?"(0). (68) 

Since it is of the form dif ~ — / 2 , we expect the solution 
R((j),£) to decay according to 1/^ and take its form to be 



(69) 



with the function Rq satisfying 

(l + io)dtRo(<f>) = R oW + \Rv{4>)Ro{& - R'^)R^(0). 

(70) 

In terms of a new flow variable t = In (I + £$) , the equa- 
tion for Ro((f>,t) is identical to that for R(4>,£), Eq. (1551) . 
with e = 1. Thus, on the scale beyond the Larkin length 



1.0 



FIG. 7: (Color online) Flow of g at d = 2: for T < vrif it 
flows to the g* = 1 — ^ fixed line, and for T > tyK, it flows 
to g* = fixed line, corresponding to a glass transition at 
T = tyK. 



(when i?o has crossed away from Gaussian fixed point 
toward the nontrivial fixed point), we find that for large 



R"(<j>,l) 



ir) 2 + 



-2 1 



18 



(71) 



with R"(0,£) = also decaying with £. This is the 
same as the result obtained by Chitra et al. [39j at the 
lower-critical dimension of d = 4 for the xy model with 
bulk random-field disorder. For this special case of an 
infinitely thick (w — > oo) 3D cell, the above result also 
reproduces the earlier finding in Ref. [2l|. 

(iii) d = 2. In two (bulk) dimensions (d — 2), O = 
leads to 7/ = T/(nK) that is fixed under the RG flow 
and the long-scale behavior is no longer controlled by a 
zero-temperature fixed point. Instead, the finite temper- 
ature selects eigenfunctions of the RG flow, Eq. (|64l) , that 
near the Gaussian fixed point are harmonics cos(r«/>) (n 
integers), with eigenvalues 



An 



1 



■ T 
ttK' 



(72) 



Focusing on the eigenfunction cos <fi with the largest 
eigenvalue, Ai, the FRG flow reduces to a standard RG 
flow equation 



dig 



1 



ttKJ 



9 - .9 



(73) 



for a single amplitude of this lowest harmonic of R(4>). 

As illustrated in Fig. [71 g decays to the g = fixed line 
for T > ttK = T g and grows to a 



1 



T 



(74) 



fixed line for T < irK = T g . Thus, we find that the 2D 
cell exhibits a Cardy-Ostlund-like [4(| phase transition 
at 



T„ = irK 



(75) 
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between a high-temperature phase, where at long scales 
surface pinning is smoothed away by thermal fluctuations 
and a low-temperature glassy pinned phase, controlled 
by the nontrivial g*(T) fixed line. This surface-pinned 
state (and the associated transition) is quite similar to 
the super-rough phase of a crystal surface grown on a 
random substrate [38j |. to the 1+1 vortex glass phase 
of flux-line vortices (confined to a plane) in type-II su- 
perconductors d, HO, El| , and to the 3D smectic liquid 
crystals pinned by a random porous environment of e.g., 
aerogel [7|. 

One key distinction here is the irrelevance in 2D of the 
tilt pinning potential [37j 



tilt 



T 



a,f3 



d d - l x{V^ - V0g) 



P\2 



(76) 



in contrast to its relevance in these other bulk pinning 
systems [3, 0, HO, HI, El| , where it leads to a super-rough 
phase characterized by In 2 x roughness. This difference 
leads to a distinct behavior of correlation functions in the 
pinned phase for the surface-pinning problem. Because 
this 2D case is somewhat academic, with our main focus 
on the experimentally relevant 3D cell, we do not explore 
it any further here. 



2. Finite-thickness cell 

We now focus on the behavior of a finite thickness cell 
with Dirichlet or Neumann boundary conditions on the 
top homogeneous substrate. Much of the RG analysis of 
the previous section extends to this case after substitu- 



tions of r^ ' and Tq V ' , Eqs. ([58)) and ((63]), respectively, 
for r(°°\ With these changes the rescalings Eqs. (|5Tj) 
and (|52[) are supplemented with the exact flow equation 
for the cell thickness 



w(b) = b 1 w. 



(77) 



The zero-temperature flow equation takes the same 
form as for an infinite cell, with a £ (T>) , (Af) , 



d t R a {$) = e^(£)R a (<b) + -R':(<f>)R'>(cf>) - R'MR'M, 



Eqs. (179a[) and (|79bl) have been, respectively, defined by 

C d -iA d ~ 3 



K 2 coth 2 (Aw) 

C d ^A d - 3 
K 2 tanh 2 (Aw) 



R(<t>), 
R(ct>). 



(80a) 
(80b) 



The limiting cases of these flow equations can be easily 
understood. For a thick cell uj —¥ oo equations for both 
(top substrate) boundary conditions reduce to an infinite 
cell analyzed in the previous section. In the opposite 
extreme of a microscopically thin cell, such that Aw <C 1 , 
the e^ a \t) functions reduce to 



c(A0 



(£) w 1 - d, for w < a, 
(£) sa 5 — d, for w <C a, 



(81a) 
(81b) 



corresponding to flow equations for a (d— l)-dimensional 
bulk random-field xy model, which in the case of the 
Dirichlet boundary condition is in a uniform external 
field. In this Dirichlet case, the eigenvalue is negative for 
any physical dimension, showing that random pinning is 
always dominated by an ordering field of the rubbed top 
substrate. In contrast, in the Neumann case, the flow as 
expected is identical to that of a (d— 1) bulk system with 
the eigenvalue of 4 — (d — 1) = 5 — d. 

For a more realistic situation of a finite cell thickness 
w, there is a crossover from an infinite (i-dimensional cell 
limit at small £ such that Aw(£) > 1 to an effective 
(d — l)-dimensional system for Aw(£) <§; 1. The cor- 
responding crossover scale is given by 6^ = w/a. An 
independent crossover scale encoded in flow equations, 
Eq. ([78]) is set by a scale b* L at which the nonlinear terms 
become comparable to the linear ones, where the flow 
leaves the vicinity of the Gaussian fixed point and (in 
the bulk system, i.e., for w — oo) would approach the 
nontrivial fixed point (|67p. From the flow equations, 
Eq. (1781 . one can see that this latter scale is simply set 
by the Larkin length, with b* L = ^l/ci- 

As we discuss below, the detailed nature of distor- 
tions strongly depends on the relative size of these two 
crossover scales and on the type of boundary condition 
on the homogeneous substrate. We naturally designate 
the two cases, w <^ £l and w 3> £l, as thin and thick 
cells, respectively. 



(78) 



except that the constant e = 3 — d is replaced by 
dependent functions e^(£), e^(£), given by 



e^(£) 



4Aw(£) 
sinh[2Aw(f) 



■, Dirichlet, (79a) 



e(A °W = e+ i lyini - Neumann, (79b) 
smh|2Aw(i:)J 

for the two boundary conditions on the top sub- 
strate. The dimensionless disorder variance functions in 



B. Correlation function 

We now turn to a calculation of correlation functions. 
As discussed earlier, on short scales (smaller than the 
Larkin length), these can be simply computed using the 
random-torque model of Sec. IIIII However, as we have 
seen m 

Sec. ELS for d < 3 the effective pinning becomes 
strong (compared to the elastic energy) on scales longer 
than the Larkin length, leading to a breakdown of the 
perturbative expansion and of the random-torque model. 
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Nevertheless, we can utilize the above FRG, which ef- 
fectively allows us to treat pinning nonperturbatively to 
overcome this difficulty. To see this, we note that the 
power of the renormalization group is that it establishes 
a connection between a correlation function at a small 
wave vector (which is impossible to calculate in pertur- 
bation theory due to the aforementioned infra-red diver- 
gences) to the same correlation function at large wave 
vectors (short scales), which can be easily calculated in 
a controlled perturbation theory 0, EH, HI, Sll • 

This relation for the Fourier transform of the surface 
(z = 0) correlation function C(q) is given by 



C[q, K, w, R a ] 



C[(ie e ,K(£),w(e),R a (£)} 



(82) 



where the prefactor on the right-hand side comes from 
the dimensional rescalings of Sec. I IV A[ remembering 
the momentum-conserving 8 function in the definition of 
C(q), and for simplicity we chose to keep T fixed under 
rescaling. We then choose the rescaling variable I* such 
that 



qe l * = A , 



(83) 



which allows us to reexpress on the right hand side of 
Eq. (|82|) in terms of the wave vector q 



C[q,K,w,R a ] = [ - 



d-l 



C[A,K(l.),w(**), ■ 
(84) 



Because the correlation function on the right-hand side is 
evaluated at the large wave vector, it is easily computed 
perturbatively in a weak pinning potential R a {4>i £*) if the 
latter is indeed small at long scale e e * a, i.e., the pinning is 
weak. To lowest order, the computation can be done with 
the replicated random-torque (Larkin) surface model 



H 



(r) 



2T(4 
which gives 

C^[A,K(Q,w(£*),R a (Q] 



£rW(4M(q)| a (85) 

OL 

1 " 

— ^K'(0,4)^(q)^(-q) 



T(4 



ri a) (4) 



ri a) (^*) 



(86) 



where in the last line we neglected the subdominant ther- 
mal part. To evaluate the resulting correlation function 



CW[q,K,W,R*] 



d-l 



ri a) (^) 



(87) 



explicitly requires an analysis of the flow for specific 
boundary conditions. 



1. Infinitely thick cell: w — > oo 



We first focus on an infinitely thick cell, for which the 
above correlation function reduces to 



C[q, K,oo,R] 



1 A d - 6 R"{0,Q 

i R"(o,e,) 

V" 1 C d . x ' 



where in the second line we used Eq. (lf?5j) to express the 
result in terms of the dimensionless disorder variance. 

(i) 2 < d < 3. As we learned in the previous sec- 
tion, for this range of dimensions, at large [which by 
Eq. (f83| corresponds to small q], the dimensionless pin- 
ning variance flows to a fixed point (|67p . giving 



C[q, K,oo,R] 



1 (3 - d)7T 2 



q 



d-l 



9C d - 



(89) 



as presented in Sec. |U 

We can now use this result to compute real-space cor- 
relations on in-plane scales x ^> £l, characterized for 
z = z' by 



C(x,z,z) = ((<Hx,z)-0(O,z)) 2 } 

= - / — ^T(l-cosq.x)e- 2 « z C7(q) 

(90) 



(27r) fl 

Cl(x, z, z) + C*(x, z, z). 



In above, Cz,(x, z, z) is a (nearly) x-independent con- 
tribution to the correlation function from short scales, 
< q < a" 1 , where Larkin approximation [random- 
torque model, Eq. (|2T)1) . from Sec.Qll] is valid and is given 
by (see Appendix IC 2\ 
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2A f f (l-cosq.x)e- 2 ^ 
/2z\ 3 ~ d 

« (3 - d)87r 2 ( — j r(d-3,2z/€i.,2z/a), for x > £ L 

- 8 7 r 2 { 1 rf ( i n - +1 - 7) ' a « 2 ^«^ fed = 2> 



C L (x,z,z) 




z 



0.5 1 1.5 2 & 

FIG. 8: (Color online) For d — 2, the contribution of short- 
scales 1/£l < q < a -1 to the correlation function C(x,z, z) 
with different a; values (from bottom to top): 0.1, 0.3, 0.5, 10, 
100£l. At small x (e.g., x — 0.1£l, blue curve at the bottom) 
Cl{x,z,z) is very small but x-dependent, while at large x 
values (shown for x — 10, 100£l, the nearly overlapping higher 
curves), the contribution is nearly x- independent: at z = 0, 
it is a constant around 8iv 2 and decays rapidly to zero. 



and is plotted in Fig. [9] 

The full 2D correlation function C(x, z, z) [denned by 
Eqs.([37|). flU])-®] is plotted in Fig. HO] 

(ii) d = 3. For a three-dimensional (infinite-thick) cell 
the pinning is marginally irrelevant, with the large I so- 
lution given by Eq. ([7Tj) . Using this inside Eq. (f88| . we 
obtain 



I 

with T(p, zi, Z2) = J z 2 t p ~ 1 e~ t dt the generalized incom- 
plete gamma function and 7 ~ 0.58 is the Euler's con- 
stant. In contrast to its small x (x <C £l) behavior, 
Eq. (|37|) , Cl(x, z, z) is x-independent for x ^> £l and is 
plotted in Fig. H 



The second long-scale part, C*(x, z,z) in Eq. (f90|) is 
a universal contribution [determined by the fixed point 
function, Eq. (|67[) ]. that for d — 2 and low temperatures 
(when 2D effects discussed in Sec. lIV Al can be neglected) 
can be straightforwardly computed. It is compactly given 
by (see Appendix IC 3|) 



(92) 



I 

The short-scale part, Cl(x, z, z), is similar to Eq. (|9"T1) . 
which when evaluated in 3D reduces to (see Appendix 
IC2I 



C*(x, z, z) 



2(3 -d)ir 2 f d d ^q 1-cosq-x 



2qz 



9C d _i J (2 7 r) d - 1 q 

1 



d-l 



9 



(2* + &) 



for x > £z,, d = 2, 
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FIG. 9: (Color online) Long-scale contribution, C„(x, z, z), to 
the 2D correlation function as a function of x,z. In the top 
(bottom) figure, it is plotted as a function of x (z) for a se- 
ries of z (x) values: z = 0, 1, 2, 10, 100£l (x = 1, 10, 100£l). 

The approximation ^ In |l + ^z+( L )' 2 ] (dotted black) pro- 
vides an excellent interpolation to the overall x, z dependence. 



Cx(x, z, z) 



2A, 



d 2 q (1 — cos q ■ x)e 



-2qz 



8tt 2 



4 Z 



r(0,2z/£ L ,2z/a), fora;>Cz 



1 



ln(2z/a) 



Stt^^ ln(5t/o) ' 

'' ' : ' : e 



a < 2z < £ L 



for <i = 3, 



(94) 



ln(a/a) 



where Cl(x, z,z), plotted in Fig. [11] is again nearly x- 
indcpcndcnt for x 3> £l, as the only x-dependence enters 
through cos q • x, which averages to zero for these large 
wave vectors with qx 1. 



For d = 3, the long-scale universal part C*(x, z, z) in 
Eq. (|90D is obtained directly from the fixed-point function 
(fTTj) [2l[ [22| , derived in Appendix IC 31 and it is given by 
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C(x,z,z) 



80 




1 2 3 4 <Tl 



FIG. 10: (Color online) Full 2D correlation function C(x, z, z) 
in the top (bottom) figure plotted as a function x (z) for 
a series of z (x) values: z = 0, 0.1, 1, 2, 10£l (from top to 
bottom) [x — 0.3,3,5, 10£l (from bottom to top)]. For z S> 
£l, the correlation is dominated by C*(x,z, z). 



C*(x,Z,z) PS —-77- 



2ir f d 2 q 1 — cosq-x _ 2qz 
9C 2 J (2tt) 2 q 2 \n(qa) 6 



2tt2 
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2vr2 
9 



1 L,„—2kz/x 



4 J ln(fca/x) 



for x > £l, d = 3, 

g— Ikz I x 



dk 



k\n(ka/x) 



k if? 7 in (2z/a) ' £l < £ < 2z 



(95) 



r 



and is plotted in Fig. [T21 

Thus on the z = substrate, at long x S> £l length 
scales and for a thick 3D cell, we find 



2tt 2 

C*(x,0,0) « —In 



ln(:c/a) 



2tt 2 

— ln[ln(x/a)], 



(96) 



as claimed in the Sec. U and first found in Ref. [2l| . The 
complete 3D correlation function C(x, z,z) is plotted in 

Fig. EH 



We conclude this section with a computation of the 
surface orientational order parameter, 



ip(w,0) = (e^( x -°)) 



-<<^( x ,0)>/2 



(97) 



where somewhat crudely we approximated it by assuming 
Gaussian correlations in (f>(x,0). This order parameter 
is of particular interest to the application of our results 
to a finite-thickness cell with a Dirichlet or a Neumann 
boundary condition imposed on the homogeneous sub- 
strate (see Fig. [1]). 
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C L (x,z,z) 




z 



0.5 1 1.5 & 

FIG. 11: (Color online) For d = 3, the contribution of short- 
scales 1/£l < q < a^ 1 to the correlation function C(x.,z,z) 
with different x values (from bottom to top): 0.1, 0.5, 1, 10, 
and 100£l. At small x (e.g., x — 0.1£l, blue curve at the bot- 
tom), (7l(x, z, z) decays very fast and is ^-dependent, while 
at large x values (shown for x = 10, 100£l, the nearly overlap- 
ping higher curves), the contribution is nearly ^-independent: 
at z = 0, it is a constant around 8tv 2 and decays rapidly to 
zero. 

Because, as we have seen above [see, e.g., Eqs. p^l) and 
(|96p]. in the limit of an infinitely thick cell (w — > oo) 
firms grows without bound with system size L, the orien- 
tational order parameter ip(w, 0) vanishes in the thermo- 
dynamic limit. For a more realistic situation of a finite 
cell, the decay of orientational order is determined by the 
cell thickness, w, and the nature of the boundary condi- 
tions on the homogeneous substrate. 

2. Crossover in a Dirichlet cell 

Contrasting the bulk behavior, for a cell with a Dirich- 
let (homogeneous) substrate, we expect the growth of 
root-mean fluctuations of <f> to be suppressed by the align- 
ment by this homogeneous substrate. Thus, in this case, 
tp(w, 0) is nonzero and the orientational order is stable 
for an arbitrarily thick (but finite) cell. 

We can analyze ip(w,0) by estimating rms (w) using 
the results of FRG found in Sec. lIV Al To this end, we ex- 
amine the asymptotics of the FRG flow in Eqs. ([77)1 and 
(1751) . For a thin cell (defined by to < £ L ), w(£) = e~ l w 
reaches the microscopic scale a at e l ™ — w/a and there- 
fore e iv \£ > t w ) « e-2 = 1- d before e 1 ^ = £ L /a. 
Since beyond t w , e^ v \£ > l* w ) < 0, pinning is irrelevant 
and the flow is cut-off at scale e £ ™ , scales beyond £l 



are not probed (the flow never leaves the vicinity of the 
Gaussian fixed point), and 4> rms can be accurately com- 
puted within the Larkin approximation (random-torque 
model), cut-off by w. 

In contrast, for a thick cell (defined by w ^> £^), the 
flow crosses over to the vicinity of the nontrivial fixed 
point i?* (leaves the Gaussian fixed point) before it is 
cut-off by the finite w. In this case, on longer scales, 
e l > w/a ^> e/ L = £z,/a, the fluctuations are cut-off 

C„(x,z,z) 




20 40 60 80 100& 



FIG. 12: (Color online) Long-scale part of the 3D corre- 
lation function C*(x, z, z) plotted in the top (bottom) fig- 
ure as a function of x (z) for different z (x) values: z = 
0, 2, 10, 100£l (from top to bottom) [x = 2, 10, 100£l (from 
bottom to top)]. The dashed curve is the approximation 
2|_ In [In (x/a)/ In (^i/a)] summarizing in in-plane correla- 
tions on the heterogeneous (z — 0) substrate. 

by w (by the flow's return to the Gaussian fixed point). 
In this thick cell regime, 4> r ms{w) is thus given by the 
matching calculation of Sec. IIVBI with the diverging L 
(> £l) dependence cutoff by w. 

Following this crossover allows us to calculate 4> rms (w). 
For a thin cell, w <C £j, and d < 3, we have 
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FIG. 13: (Color online) Full 3D correlation function C(x, z, z) 
in the top (bottom) figure is plotted as a function of x (z) 
for a series of z (x) values: z = 0,0.01,0.2, 1,2£l (from top 
to bottom) [x = 0.3, 3, 5, 10£l (from bottom to top)]. For 
z 2> £l, the correlation is dominated by C**(x, a, z). 



jd-1, 



" 1 J^y 

w ' 1 d^q A f w 2 d d - x q A f 



(2w) d - 1 K 2 J w -i {2n) d - 1 K 2 q 2 

3-d 



« 4tt^ I —J , forw«a,rf<3, (98) 

where we used the Larkin approximation Eq. (15"3")) together with Eq. , valid for w <C £j, since the Dirichlet flow in 
Eq. |78|) never leaves the vicinity of the Gaussian fixed point. 
For a thick cell w ^> £l and d < 3, we have 

jci-l 



" 1 n VD 



(27T) 

^ d-*-^ 67T 2 1 r' 1 d d - x q A } 



-1 (2ir) d - 1 9C d -! q d - x J^-i (2ir) d - 1 K 2 q 2 
eir 2 

— ln( W /a) + 47r 2 , for w > rf < 3, (99) 
I 



where we neglected the subdominant contribution from around" heading back toward the Gaussian fixed, i.e., 
scales longer than w (where the RG flow for R(£) "turns pinning is irrelevant), in the first term of the second line 
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approximated the <j) correlator C{q) by its fixed-point 
value valid for w/£l >• 1 (such that the flow ap- 

proaches the vicinity of the nontrivial fixed point), and 
approximated C(q) in the second term by its Gaussian 
fixed point expression (Larkin approximation), valid for 
£jT < q < a -1 . We also used the definition of £l to 
approximate the second term in the last line by (27r) 2 . 

Repeating above estimates for d — 3, for a thin cell 
(w <C we find 



4-7T 



d 2 q A 
i (2^)2 lOq 2 
2 ]n(w/a^ 



for w<c£ L , d = 3, (100) 



be enhanced compared to the w — > oo system, where ad- 
ditional homogeneous bulk degrees of freedom have a sta- 
bilizing effect against pinning (c£/ c reduced from 4 down 
to 3). 

This is reflected in the behavior of both the correlators 
in the random-torque model [given by Eqs. (f2"5"j) and (155")) ] , 
and in the FRG flow that becomes only more divergent 
on scales e l > e l ™ = w/a, as e {N \l > £* ) ->• e + 2 = 
5 — d > e. Hence, in contrast to the Dirichlet cell (where 
finite w suppresses the effect of the random potential), 
in a Neumann cell, finite thickness enhances the effects 
of random surface pinning. Consequently, independent 
of the Neumann cell thickness, the orientational order 
parameter, ip, vanishes for L — » oo. 



and for a thick cell (w 3> 



STRONG PINNING LIMIT 



d 2 q 7T 2 



-i (2tt) 2 9C 2 q 2 ln(qa) 



d 2 q A f 
J^flOq~ 2 



hi 



ln(w/a) 
ln(a/a) 



47T 2 , for w > £ L , d = 3, 



(101) 



where we employed the same asymptotic approximations 
as for d < 3. 

Putting these crossovers together, we finally obtain the 
surface orientational order parameter ip(w, 0) for thin 
and thick cells in d < 3, 




thin cell, w -C £l 



" 2 thick cell, w> a, 



(102) 



thin cell, w <C £l 
thick cell, w >• £l, 



where 77^ = (3 — g?)7t 2 /18, 77313 = 7r 2 /18 are univer- 
sal exponents [given in Eqs. (|99| and dlOip ] and i]l — 
27r 2 /ln(^/a) is a nonuniversal constant. The 3D sur- 
face order parameter for such Dirichlet cell of thickness 
w is illustrated in Fig. O 



3. Crossover in a Neumann cell 

Above analysis straightforwardly extends to a finite- 
thickness cell with the Neumann boundary condition on 
the homogeneous substrate. At long scales, a finite- 
thickness Neumann cell reduces to an effective "film," i.e., 
ad— 1-dimensional bulk random-field xy model. Thus, 
we expect the disordering effect of the random pinning to 



In all of the above analysis, we focused on the most 
interesting weak surface disorder, where pinning is col- 
lective, dominating over the elastic energy only on the 
macroscopic length scales, longer than £l 3> a. This 
assumption is what justified our treatment of the elastic 
energy as dominant (at least on short scales, smaller than 
£l), allowing an expansion about the ordered </3 = (ne- 
matic) state. However, it is quite possible that in some 
(e.g., liquid crystal) applications, it is the opposite limit 
of strong pinning that is of interest. 

In the latter strong-disorder limit, the surface-pinning 
potential (by definition) dominates over the elastic en- 
ergy at all, even microscopic scales, with V p > K/a. To 
treat this regime, we instead perturb in the elastic energy 
about a random ground state, </3g(x), that exactly min- 
imizes the random pinning potential V[0o( x )j x ]- That 
is, (90 o V r )|^ o ( x ) = 0s( x ) = 0. We then expand about this 
random ground state, obtaining 



jd-1 



9 



(2tt 



fg|^o(q)| 2 + f|0o(q)-^(q) 



(104) 



(103) wh ere g = -d 2 o V\ M ^ )=4> s {x) « V p /(2tt) 2 . A minimiza- 



tion of the above Hamiltonian then straightforwardly 
gives 

.9 



0o (q) 



Kq 



:^o(q)- 



(105) 



From this analysis we can readily identifying a strong- 
coupling pinning length 



K/g, 
K 



{d-l)/2 



A 1/2 ° 

A3-d)/2Ad-l)/2 
SL SO ' 



(106) 



which is a scale below which the xy-order parameter no 
longer faithfully follows spatial variations of the local ran- 
dom potential and thus it is a strong-coupling version of 
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the Larkin length. In above, we restored the pinning 
potential correlation length £ to also account for the 
more realistic case where £o is distinct and longer than 
the microscopic molecular cutoff scale a. Hence we con- 
clude that below the pinning correlation length £o , there 
is a crossover from weak to strong pinning limit when 
the (weak-coupling) Larkin length, £l, drops down to £ s . 
On these shorter scales, the collective pinning analysis of 
previous sections and corresponding results break down. 



VI. APPLICATION TO LIQUID CRYSTAL 
CELLS 



As discussed in the Sec. HI liquid crystal cells provided 
a strong motivation for our study of the orientational 
order in the presence of surface random pinning. How- 
ever, although there is a qualitative overlap, in detail a 
model of a surface-pinned liquid crystal cell, a priori can 
be quite different from a basic xy model studied above. 
Furthermore, the detailed model very much depends on 
the specific nature of the liquid crystal phase and thus 
requires an extensive study that lies beyond the current 
paper. However, to put our above results for an xy model 
in a physical context, we now briefly examine a formula- 
tion of a surface-pinning problem for real liquid crystals, 
focusing on nematic and smectic phases, deferring their 
detailed analysis to a future study. 



A. Nematic liquid crystal phase 

The key distinction of the nematic liquid crystal phase 
as compared to the xy model studied so far is the nature 
of the Goldstone modes, that for a nematic is given by 
a i/iree-dimensional unit vector (strictly speaking with 
opposite ends identified forming an RP2 manifold), the 
nematic director 



(cos 9 cos cp, cos 9 sin ip, sin ( 



(107) 



as opposed to its xy-model counterpart, where it is a sin- 
gle azimuthal (planar) angle, (j>, with the polar angle 9 
fixed at zero by some easy-plane anisotropy. We note 
that in above, we chose a somewhat nonstandard (but 
here convenient for treating parallel surface alignment) 
convention for 9. We also implicitly ignored the diffi- 
cult question of topological defects proliferation and the 
corresponding stability of the elastic glass. If indeed im- 
portant, for weak disorder, we expect their effects to set 
in on much longer length scales, thereby providing a wide 
intermediate range of scales, where defect-free model is 
of interest. 

In the nematic phase, the bulk energy of a nematic 
director n(r) is described by the well-known Frank-Oseen 
expression [43| 



H F = i J d d - l xdziyK s (V -h) 2 + K t [h-{S7 x n)f 



+K b [n x (V x h)] 



(108) 



Within a simplifying one-elastic-constant approximation, 
K s = K t = Kb = K n , and together with the surface- 
pinning energy and polar parameterization, Eq. (11071) . 
the above equation gives the elastic Hamiltonian for the 
nematic surface-pinned cell 



tt _ f^n I id- I 

-^nematic 2 



J d d - x x J dz [(V0) 2 + cos 2 9{Vip) 2 } 
+H pin , (109) 
where the surface-pinning energy 



TT I jd-l. 



+W W COS 



(W Q + V(2v,x))cos 2 9\ z=0 



W v I Z — W 



(110) 



is given by a purely homogeneous planar (9 = 0) align- 
ment on the top (z = w) substrate and a planar align- 
ment with a heterogeneous azimuthal component on the 
bottom (z = 0) substrate. The latter is encoded in a ran- 
dom pinning function V((j>, x) with a 2ir periodicity of </>, 
presented in Sec.|Hl with cj> = 2ip capturing the h <-> — n 
symmetry of the nematic liquid crystal phase. 

For a thin cell and weak random pinning on the bottom 
substrate, such that the scale of V((f>, x) is much smaller 
than Wo, clearly the planar alignment, while of random 
strength (on the bottom substrate), remains planar on 
both substrates. We further note that because (V</>) 2 
(computed within the xy-model approximation) remains 
finite in the physically interesting dimensions (i.e., in- 
frared convergent for d > 1) and small for weak pin- 
ning (decaying into the bulk with z), based on Fredericks 
transition phenomenology |43| . we do not expect planar 
surface alignment by Wq. w to be overturned by a weak 
second term in Eq. (|109l) . 

Thus, with the exception of small regions (that are rare 
for weak disorder), 9 — is the solution that minimizes 
the total energy. With this, the nematic cell model re- 
duces to the surface random-field xy model for the di- 
rector's azimuthal orientation 4> = 2(p, studied above 
[2H . l22j . Hence, for weak pinning, all of a;y-model re- 
sults detailed above apply directly to a nematic liquid 
crystal cell. 

In contrast, we expect strong pinning to lead to large 
azimuthal distortions that will be accompanied by big 
9 variations both on the random substrate and in the 
bulk. For intermediate pinning strength, the system can 
perhaps even exhibit a random Fredericks-like transition 
corresponding to a bulk escape from a planar surface con- 
figuration, driven by a large (V0) 2 3> 1/w 2 . We leave 
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the detailed study of the associated subtleties for this 
system, which distinguish it from the simple xy model to 
a future research. 



B. Smectic liquid crystal phase 

Another important realization of a random surface- 
pinning problem is that of a smectic liquid crystal on 
a heterogeneous substrate, as for example realized in re- 
cent experiments [26| mentioned in the Sec. |H We will 
focus on the experimentally and theoretically more in- 
teresting case of the bookshelf geometry, illustrated in 
Fig. [TJ where layers and director are, respectively, per- 
pendicular and parallel to the substrate. Choosing the 
coordinate system as indicated in Fig. HH so that the 
smectic layers lie parallel to the (x, y) plane, with the 
average layer normal along the z axis and the random 
substrate located at y = and running perpendicular 
the y axis, the total energy is given by 



Ha 



= a 



dy 



§(^) 2 



+H p 



(111) 



where the first two terms describe the usual smectic elas- 
ticity, with bending and compressional elastic constants 
K and B (for simplicity taken to be harmonic [3), and 
the last term is a surface-pinning energy given by 



d d - l xdy5{y) 
d^xdySiy) 
d^xdySiy) 



W 



(h ■ y) 2 - (n ■ g(r)) 2 - V[u,r) 



-—(Sn y ) - h(x)5n x - V[u,x] 



W 

— {dyii) - h(x)d x u - V[u, x] 



(112) 



where for convenience we defined r = (z,x,y) = (x, y) 
and extended its x to d — 1 dimensions transverse to y. 

In above, the W and g(x) terms are the homogeneous 
and random components of the orientational pinning [W 
inducing a homogeneous parallel to the surface director 
alignment and g(x) capturing random azimuthal director 
pinning within the heterogeneous substrate plane] and 
V[u, x] is the positional pinning of surface layers with 
the random substrate [1, 0| • In getting to the final form 
we expanded h about its pinning-free orientation along z 
and used the smectic "Higgs mechanism" 0, Hl| to make 
a replacement dn± — > V^u, valid inside a smectic phase. 

Following procedure used for the xy model, this ran- 
domly surface-pinned smectic model can be used to an- 
alyze the phenomenology of a smectic liquid crystal cell 
with the hopes of understanding long-scale random tex- 
tures observed in Ref. [26[ . Because (as illustrated above) 
the smectic elasticity and pinning differ qualitatively 
from that of an xy model, we expect a phenomenology 




FIG. 14: (Color online) Cartoon of smectic liquid crys- 
tal cell with random substrate at y — 0. The substrate 
fixes large orientation transversely, modeled by W(Sn y ) 2 ~ 
W(d y u) 2 , and the surface disorder randomly pins layer po- 
sitions through V[u, x]5(y) and layer orientations through 
(n • g(r)) 2 ~ h(x)d x u. 



that is qualitatively distinct from that found for the sim- 
ple xy model and the nematic phase found above. We 
leave the interesting and nontrivial study of a random 
smectic cell to a future research. 



C. Experimental observables 

One attractive feature of liquid crystals is that their 
orientational order can be readily studied via light mi- 
croscopy. In its simplest form, the technique utilizes a 
crossed polarizer-analyzer pair on the front and back of 
a cell, typically transversely oriented. In this geometry, 
the spatial (within the xy plane) distribution of the trans- 
mitted light intensity through the cell is sensitive to the 
azimuthal variation of the local optic axis and therefore 
measures the director's planar spatial distribution. 

For a fully ordered planar nematic state, with a uni- 
form director orientation at an azimuthal angle <p with 
respect to the polarizer (or analyzer), the transmitted 
light intensity / through a uniaxial cell of thickness w is 
given by 43] 



I w = /osin 2 (2(p)sin 2 0V2). 



(113) 



In above, S = \o — Xe = 27r(n — n e )w/X is the phase dif- 
ference between ordinary and extraordinary components 
of light with wavelength A, respectively, characterized by 
n Q , n e indices of refraction. In the simplest case of the 
director uniformly aligned along the polarizer (ip — 0) or 
the analyzer (ip — 7r/2), this leads to a uniformly van- 
ishing transmitted light intensity. Conversely, the max- 
imum light transmission is produced for a 7r/4 uniform 
director orientation relative to the polarizer (or equiv- 
alent^, transversely-crossed analyzer). The last factor 
in Eq. (II 13[) leads to transmission color selectivity with 
optical anisotropy and cell thickness. 
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For a spatially nonuniform director variation, the anal- 
ysis of the transmitted light intensity is more compli- 
cated. However, for a nematic variation on a scale longer 
than light's transverse coherent length (a typical situa- 
tion for illumination with incoherent source), the output 
intensity simply images the transverse (xy-) optic axis 
variation, with each coherence region treatable as an in- 
dependent column of depth w. Furthermore, in the limit 
that spatial variation along z is also smooth on the scale 
of light's wavelength, the transmission through each col- 
umn can be treated in Mauguin limit, where light com- 
ponents along and perpendicular to the optic axes simply 
adiabatically follow the local director orientation n(x, z). 

For a 7r/2-crossed polarizer-analyzer pair, standard 
analysis in this Mauguin limit then gives the output light 
intensity (behind the analyzer, at z = w) 

I w (x) = Iq\ cos(^(x, 0) sin<^(x, w)e lXc 
— sin ip(x, 0) cos <p(x, w)e lXo | 2 
= I sm 2 [ip{x,w) - v3(x,0)]cos 2 (<V2) 

+/osin 2 [^(x,w) +^(x,0)]sin 2 (5/2). (114) 

For the director on the back substrate aligned (by the 
Dirichlet boundary conditions) with the polarizer axis 
(and perpendicular to the analyzer), i.e., (p(x,w) = 0, 
the output signal simplifies considerably to 



I w (x) = 7 sin 2 [^(x,0)] 



(115) 



and is thus directly related to the local surface orien- 
tational order parameter, ?/>(x, 0) = e* 2 ^*' ) = e l( ^ x '°) 
studied in this paper. For example, a spatially averaged 
transmission through a Dirichlet cell is given by 



(Ux)> 

h (l - <e**>(*))) 
~/ (1-^(«;,Q)), 



(116) 



in which ip(w,0) is computed in Eqs. (|102l) and (|103l) . 
Thus, for this choice of geometry, a thin Dirichlet cell has 
the expected vanishing transmission, that grows with cell 
thickness to its maximum value of 1/2. More stringent 
tests of our predictions can further be made by com- 
paring (xy-) spatial correlations of light transmission, 
(I w (x)I w (x')), with orientational correlation functions to 
which these are clearly directly related according to 



(I w (x)I w (x>)) 



7 2 (sin 2 ^(x,0)sin 2 ^(x',0)) 
'1 1 



-Wg>/2 



+ l e -((0o(x)-0o(x')) 2 >/2 



(117) 



An even more direct probe of director correlations is 
possible through the polarized confocal microscopy [45| . 



where an image of the local director orientation at each 
depth z can be produced. A numerical computation of 
thereby measured director correlation functions therefore 
allows a detailed comparison to results predicted here. 



VII. SUMMARY AND CONCLUSION 

In this paper, we have studied the stability of random 
distortions in an xy model perturbed by a random sur- 
face pinning and discussed our findings in the context of 
nematic liquid crystal cell with a dirty non-rubbed sub- 
strate. We found that for a thick 3D cell, at long scales, 
the disordering effects of the random substrate always 
marginally dominate over the bulk nematic order. Thus, 
a 3D nematic order is marginally unstable with orienta- 
tional "roughness" growing as ln[ln(x/a)] on long scales. 
We have also extended these results to a finite-thickness 
cell, with a second homogeneous substrate with parallel 
Dirichlet and Neumann boundary conditions. Not sur- 
prisingly, in the former case, the nematic order is stabi- 
lized to arbitrary long scales, but with the nematic order 
parameter (and the corresponding birefringence) exhibit- 
ing a crossover from a large value for a thin (weakly het- 
erogeneous) cell to a small value for a thick (strongly 
heterogeneous) cell at a characteristic cell thickness set 
by the Larkin length, £j> We expect our predictions to 
be experimentally testable via a polarizer-analyzer trans- 
mission microscopy and by studying how the nematic or- 
der is recovered in response to a tunable in-plane aligning 
electric or magnetic field. We propose that the predicted 
statistical properties (correlation functions) of the ran- 
dom substrate-induced director textures can be quanti- 
tatively tested with the polarized confocal microscopy 
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Appendix A: Larkin lengths analysis 

In this appendix, we provide the details for the anal- 
ysis of the Larkin length [9| in finite thickness, two- and 
three-dimensional cells. As derived in Sec. lIIIl the Larkin 
length is defined in the standard way, given by 
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Finite thickness in two dimensions 



6§(x)> = (27T) 2 = 
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qA f 
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with the r^ a) 's given by Eqs. ([2"3" |) -P5 |) and the lower 
momentum cutoff of the integration given by 1/£l- in 
the limit of an infinitely thick (w = oo) cell above integral 
is straightforwardly computed, in 2D (d = 2) giving 



A/ r dq 

- *k? J mr q 2 

A/ 



too 

which leads to a Larkin length 
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with the superscript oo denoting the result of an infinitely 
thick cell (that for simplicity of notation we will drop). 

In 3D, with the two-dimensional random substrate, we 
obtain 



<$(X)) = (27T) 2 
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d 2 q 1 
K 2 J (2tt) 2 'q 2 
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which gives the Larkin length as 

<- v / AT 

S.L.3D = ae f . 

For a general dimension d < 3, we have 
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giving 
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where we ignored the strongly subdominant (for d < 3) 
a term. 



For a finite-thickness (u;) 2D Dirichlet cell, the surface 

(T>) 

variance determines the Larkin length t; L according to 

''" (A8) 
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By scaling variables, this defines an implicit expression 

1 



(If')- 1 y 2 coth 2 (y) 



dy = £,l, 



for the Dirichlet Larkin length S,i°\w,^l) — 
to^ (£,l/w) in terms of the infinite cell's Larkin 
scale £l, latter simply a characterization of disorder 
given by Eq. (|A"3i 

Similarly for a 2D Neumann cell, we have 
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which gives 



1 



-dy = Z L , (All) 



i(£M)-i y 2 tanh 2 (y) 
both evaluated numerically and plotted in Fig. [51 

2. Finite thickness in three dimensions 

Repeating the analysis in 3D for the Dirichlet cell, with 
scaled ultraviolet cutoff a = a/w and 
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gives 
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which reduces to 
»i 



dy 



-= In (1.18a). (A14) 
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To get to this final result, we used the fact that the second 
integral in Eq. (|A13I) is finite in the ultraviolet, for w S> a 
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giving an a- independent constant about —0.17, thereby 
eliminating dependence on a. 

For the Neumann cell, we instead find 



with 



dy 



)-i ytanh y 



ln(0.79£ L ), (A15) 



where again a physically relevant limit w 3> a was taken 
to eliminate the a dependence. Numerical evaluation 
of Eqs. (|A14|) and (|A15[) gives the 3D results shown in 
Fig. E 



3. Larkin length crossover 



As we can see from Eqs. (|A"9l) . (lAllj) . (IA14P and (|A15I) . 
the Larkin length in a cell of thickness w depends on a sin- 
gle dimensionless ratio, of the infinite cell Larkin 
length (characterizing pinning strength) to the cell thick- 
ness. We expect that for a thick cell (w ^> the re- 
sult of infinite thick cell should be recovered. On the 
other hand, for thin cells (w <C we expect the ho- 
mogeneous substrate boundary condition to play a role. 
Namely, since the Dirichlet boundary condition on the 
top homogeneous substrate explicitly orders the director, 
suppressing the distortions of 4>, we expect £l to diverge 
for a thin Dirichlet cell. Furthermore, since the Neumann 
boundary condition eliminates the stiffening by the bulk, 
in the thin Neumann cell we expect £l to approach the 
value for a (d — l)-dimensional bulk system with (d — 1)- 
dimensional pinning. The expected crossover is indeed 
confirmed by a numerical evaluation with the solution 
illustrated in Fig. [BJ with the Dirichlet £l diverging at 
£l/w « 1.71 in 2D and £ L /w m 1.23 in 3D. 

The divergent asymptotic behavior can be obtained by 
expanding the implicit expression for 0, in Eqs. (|A9|) and 
(|A14|) around 1/£ L = 0. For d = 2, we have 
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which making use of —* — A 
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with 
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For d = 3, using ycot ^ (y) -> y as y -> 0, we have 
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& ~ 1.233. 



This leads to 
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quoted in the main text and consistent with the numeri- 
cal evaluation of the integral solution. 



Appendix B: The cusp and the fixed point 

For completeness, we now fill in some of the details 
(previously reviewed in Ref. (46|) for the RG evolution 
of the random potential variance R{4>) into its universal 
cusped form. To this end, by differentiating the flow 
equation (po) for R((f>, £) with respect to <j), we obtain 

d e R"(<f>) = eR"{^) + i?"'(0) 2 + R"(cj>)R""{<f>) 

-R""(<j>)R"(0), 
dtR""(4>) = e^""(0)+3i?""(0) 2 +4 J R"'(0)i?( 5 )(0) 

+R"{0)R^((f>) - i? (6) (0)i?"(O). (Bl) 

Setting <f> to 0, we obtain 

d e R"(Q) = ei?"(0) + i?"'(0) 2 ^ei?"(0), 
diR""(0) = eR""(0) + 3R""(0) 2 + 4i?"'(0)i? (5) (0) 

-> eR""(0) + 3R""(0) 2 . (B2) 

We note that R(<j>) is an even function and moreover 
(before the cusp develops) is smooth in </>, with R"'(0) = 
fl(5)(0) = 0. 

Clearly, the flows of R"(0) and R""(0) are exact, with 
the later diverging in a finite RG time according to |ll| 
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1 - 3c(e rf -l)/e' 



(B3) 



with c = R"" (0)|^ = o. For a special case e = (d = 3 
in our system with surface pinning), the cusp develops 
according to 



R""(0)U 



l-3c£ 



(B4) 



Thus after a finite RG time £, R""(0) diverges, signal- 
ing the appearance of a cusp in R" (</>), as illustrated in 

Fig. rrsi 
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Renormalization 



-R'X<P) 



-R"(</>) 




3n-2n-n n In 3n -3n-2n-n n 2n 3n 



FIG. 15: (Color online) Evolution of R"(<j>) and appearance 
of cusps under coarse-graining renormalization. 



Appendix C: Details of asymptotics of various 
correlation functions 



As derived in the main text, in an infinitely thick 
cell, the momentum space correlation function has short- 
and long-scale limits. The former one, computed in 
the random-torque (Larkin) approximation (scale shorter 
than Cl) on the random substrate (z = 0), is given by 



C{q) 



A 



K 2 q : 



, q > VCl- 



(Cl) 



In the long-scale limit (small q < 1/Cl), derived via FRG 
and matching methods in Sec lIV Bl it is instead given by 



9 q z \l\(qa) ' ' 



A Fourier transform of C (q) [using above limits and gen- 
eralized to finite z; see Eqs. (p~8|) - (f20|) ] then gives this 
correlation function in real space. 



1. Mean-squared distortion in Larkin 
approximation 

The mean-squared distortion of </>(x, z) can be given as 

jd-i. 



<^(x,z)>* ; / , d t f , e~^ 



(27r) d - 1 K 2 q 2 
For d = 2, we have (with q = 2qz) 
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in which the integral converges at large q so we can ignore 
the upper cutoff A = \/a of q and the definition of Cl as 
in Eq. ^ was used. Expanding this result at different 
range of z values, we obtain 

' 1 - |f( ln 27 + 1 ~ 7), a < 2z < £l 



I 2z 
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in which 7 « 0.58 is the Euler's constant. 
For 2 < d < 3, we have 
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making use of q = 2qz and the definition of £t as in 
Eq. (J3J. For different range of z values, we obtain 
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Similarly, for d = 3, we need to consider the upper 
cutoff of q and have 
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where g = 2gz. The mean-squared distortions for d = 2 
and d = 3 are plotted in Fig. [5] 



2. Correlation function in Larkin approximation 

The short scales (where Larkin approximation holds) 
contribution to the real-space correlation is given by 
Eq. (|36|) . In an infinite thick cell, we have 
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in which the integral of q has lower cutoff 1 /Cl and upper 
cutoff 1/a. For simplicity of notation, we will ignore the 
superscript 00. 

For d = 2 and x <C Cl, since the kernel is convergent at 
large q, we can extend the integral to infinity (requiring 
x 3> a) and obtain 
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On the heterogeneous substrate (z — 0), we have 
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where Si(z) — Jg s " 1 t ( - t - ) dt is the sine integral function. 
Making use of the expansion of cos (qx) at qx <C 1, for 
finite z we obtain 
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Stt 2 / f 1 / 
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where the last approximation is taken for z»i and the 
other terms are subdominant. 

When x ^> £l, the cos(qx) oscillates strongly giving 
subdominant contribution. Then we have 
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in which the generalized incomplete gamma function 
T(p, zi,z 2 ) = t p ~ 1 e~ t dt — T(p, zx) — T(p, z 2 ) and we 
made use of expansions T(— 1, x) ~ ^ + (7— 1+lnx) — x/2 
at small 1 (1 < 1) and r(— l,x) ~ x~ 2 e~ x at large x 
(x> 1). 

For 2 < d < 3 and x <C £l, we can extend the integral 
to infinity (requiring x 3> a) and obtain 
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in which fd{qx) is the average of cos (q • x) over a sur- 
face of a (d — l)-dimension unit sphere. The correlation 
function could be evaluated numerically and on the het- 
erogeneous substrate this result approaches (with q = qx) 



C L (x,0,0) 



8(3 -dK 2 r i-fM 
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as given in Sec. IIII1 Making use of the expansion 
1 — fd{qx) ~ 2(d-i) ^ ^ ^ 1, we can evaluate the 



correlation function at finite z as 
C L (x, z, z) 
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in which the last approximation is taken for z 3> x and 
the other subdominant terms are ignored. 

When x 3> £l , the cos (q ■ x) again oscillates at the 
full range of q and obtain 
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where we made use of expansions r(d — 3, x) ~ T(d — 3) + 
at small x (x <C 1) and T(d — 3,x) s» x d_4 e _;l: at 
large x (x ^> 1). 

For d = 3, the variables q and x are two dimensional, 
thus we have 
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in which Jo(qx) is the Bessel function of the first kind. 
The behavior of its correlation could be evaluated nu- 
merically and the asymptotics can be obtained approxi- 
mately for different regions. 

On the heterogeneous surface (z = 0) , for small x (a <C 
x <C we have 1 — Jo(qx) s» 1 with q ^ l/x and 
1 — Jo(qx) « (qx) 2 /4 with q -C l/x, so in this region we 
have 
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Making use of the expansion of Jo(qx), we can evaluate 
the correlation function at finite z as 
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where in the last approximation we kept the leading term 
for z^> x. 

When x~^> the expression simplifies to 
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in which we made use of expansions r(0, x) ps (—7 — 
In a;) + a; at small x (x -C 1) and T(0,x) ps x _1 e _a: at 
large a; (a; >• 1). 



3. Universal (long-scales) part of correlation 
function 



The correlation function in momentum space at small 
q (q < is obtained by FRG and matching methods 

in Sec. IIVB1 Here we calculate the corresponding real- 
space correlation functions. By construction, this form 
of C(q) only holds at < q < , with f£ therefore 
entering as the upper (UV) cutoff on all q integrals done 
here. 



For d < 3, the FRG derived correlation function rated 
— ] 
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at q < £ r 1 is given by 
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A better approximation is obtained by using a "soft" 
upper cutoff by inserting a factor of e~ q ^ L inside above 
integrand. 



For d = 2, this correlation function is given by 
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where L is system size to be taken to 00 at the end 
of calculation. Using the gamma function expansion 
r(0, x) ps (—7 — In a;) + ar, we obtain 
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Except for approximations associating with the match- 
ing method, this real-space result is an excellent approx- 
imation to a numerical integration of C(g), as shown in 
Fig. [SJ On the heterogeneous substrate and for x ^> £l, 
it reduces to C*(x, z, z) pa 2|_ m (xj^ L ). 

For <i = 3, the correlation function is given by 
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which we evaluated numerically with a soft cutoff 
i+(q£l)' 2 anc ^ Pitted in Fig. [T3] 

Approximate asymptotic behavior of this correlation 
function can be obtained analytically 
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in which k — qx. 

On the heterogeneous substrate (z = 0), the first inte- 
gral is negligible, thus for x ^> £l we have 
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Further away from the substrate such that £l <C 2z -C 
x, the e _29Z acts like upper cutoff at q ~ l/2z, giving 
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Finally, for <C x <C 2z, we find 
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with numerical prefactor that is off by a factor of i rela- 
tive to the numerical integration that gives Cj|s (x, 2, z) ~ 

? ln(2 1 Z /a) ' aS § iven m SeC - 
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